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ABSTRACT 

We show the effectiveness of strong lensing in the characterisation of Lyman continuum 
emission from faint L ^ 0.1L* star- forming galaxies at redshift £ 3. Past observations 
of L L* galaxies at redshift £ 3 have provided upper limits of the average escape 
fraction of ionising radiation of f esc ~ 5%. Galaxies with relatively high f esc (> 10%) 
seem to be particularly rare at these luminosities, there is therefore the need to ex- 
plore fainter limits. Before the advent of giant ground based telescopes, one viable way 
to probe f esc down to 0.05 — 0.15L* is to exploit strong lensing magnification. This 
is investigated with Monte Carlo simulations that take into account the current ob- 
servational capabilities. Adopting a lensing cross-section of 10 arcmin 2 within which 
the magnification is higher than 1 (achievable with about 4-5 galaxy clusters), with a 
U-band survey depth of 30(30.5) (AB, la), it is possible to constrain f esc for z ~ 3 
star-forming galaxies down to 15(10)% at 3-sigma for L < 0.15L* luminosities. This 
is particularly interesting if f esc increases at fainter luminosities, as predicted from 
various HI reionization scenarios and radiation transfer modelling. Ongoing observa- 
tional programs on galaxy clusters are discussed and offer positive prospects for the 
future, even though from space the HST/WFC3 instrument represents the only op- 
tion we have to investigate details of the spatial distribution of the Lyman continuum 
emission arising from z ~ 2 — 4 galaxies. 

Key words: galaxies: distances and redshifts - galaxies: high-redshift - gravitational 
lensing: strong 



1 INTRODUCTION 

Recent works suggest that at z > 3 star-forming galaxies 
are the leading candidates for the production of ionising 
photons (e.g., Kuhlen & Faucher-Giguere 2012; Haardt & 
Madau 2011; Ciardi et al. 2012). However, the mechanisms 
regulating the escape fraction of ionising radiation (/ esc ) 
from galaxies are still unknown. In particular it is not clear 
if / esc evolves with cosmic time, and at fixed redshift, if 
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it is luminosity dependent. The complexity in modelling of 
galaxy evolution and the inclusion of radiative transfer pre- 
scriptions make the predictions on / esc very uncertain, and 
opposite results are often obtained from simulations. More- 
over, the predictions on / esc show a large variation from 
galaxy to galaxy, between 0.01 to nearly 1, as a result of 
differences in the hydrogen distribution (e.g., Gnedin et al. 
2008; Yajima et al. (2011); Fernandez & Shull (2011), and 
references therein). From an observational point of view, 
the current situation is also far from clear. Various obser- 
vations have provided significant upper limits on Lyman 
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continuum (LyC) emission in the redshift range 1-4, i.e., 
f esc smaller than 5-10% (Vanzella et al. 2012, V12, and ref- 
erences therein). In particular Vanzella et al. 2010, VlOb, 
and Boutsia et al. (2011) reported upper limits on f esc of 
5% for LfiL* LBGs at redshift 3-4. Very few galaxies with 
possible LyC detection have been reported in the literature 
(Shapley et al. 2006; Iwata & Inoue 2009; Nestor et al. 2011). 
The fact that bright sources with high f e3C (> 10%) are 
rare could support the interpretation that f e3C increases at 
lower luminosities. Indeed, recent radiative transfer calcula- 
tions coupled with cosmological simulations show that low 
luminosity and dwarf galaxies have high f esc (> 20 — 40%) 
and are the major contributors to the ionizing background 
at redshift 3 < z < 6 (e.g., Wise & Cen 2009; Razoumov & 
Sommer-Larsen 2010; Yajima et al. 2011; Fernandez & Shull 
2011; but see Gnedin et al. (2008) for different results). The 
investigation of faint (and/or low mass) star- forming galax- 
ies at moderate redshift (z ~ 3 — 4) is therefore crucial to 
probe ionisation regimes and conditions that would be in 
place during the reionization epoch (z > 7). It is worth not- 
ing that the reionization at z = 7 — 10 requires an average 
/ esc much higher than what is observed for relatively bright 
and lower redshift galaxies (z < 4) (Robertson et al. 2010; 
Bouwens et al. 2010; Fernandez & Shull 2011). 

While the direct measurements of escaping LyC photons 
are prohibitive during the epoch of reionization, they are still 
accessible at z ,$ 4. The investigation of sub-L* galaxies re- 
quires extremely deep U-band surveys (U > 30). However, 
even if they were available, spectroscopic redshift measure- 
ments for L « L* galaxies (e.g., r-band> 27) are very 
challenging with current facilities. In this respect, before 
the advent of the ELT and JWST telescopes, strong lensing 
magnification offers a viable way to explore faint limits in 
luminosity (and mass). 



2 THE f esc AND STRONG LENSING 

The relative fraction of escaping LyC photons (at A < 912 A) 
relative to the fraction of escaping non-ionising ultraviolet 
(1500 A) photons is defined as (Steidel et al. 2001): 

/csc - rcl = Te exp(r LHC ), (1) 

V F LyC > ohs 

where (F1500 /Fl v c) obs, (£1500/ LL y c)int and represent 
the observed 1500A/LyC flux density ratio, the intrinsic 
1500A/LyC luminosity density ratio, and the line-of-sight 
opacity of the IGM for LyC photons, respectively. It is worth 
noting that the two quantities (fT500) o bs and (FL y c)obs are 
measured in the same spatial (i.e., physical) region, where 
the ionising and non-ionising radiation arise (V12). If the 
dust attenuation ^1500 is known, / eS c,rei can be converted to 
the absolute / csc as / csc = 10~°' 4Al50 ° / CS c,rci (e.g., Siana et 
al. 2007). The term /1 is the magnification factor provided 
by a lens (e.g., a cluster of galaxies) and applies to both 
flux densities observed at two wavelengths (we leave the pa- 
rameter explicitly in the expression). We assume for sim- 
plicity it does not vary substantially within the isophote of 
the sources. Therefore, given the achromatic nature of lens- 
ing, f e3C is independent from fi. The two quantities / esc and 
fesc,rei are equal if ^1500 = 0. Conservatively, in the follow- 
ing we perform calculations assuming ^1500 = 0, the limits 



probed on f esc are deeper if A1500 > 0, e.g., they are halved 
if A1500 = 0.6 (VlOb). If not specified, we assume an in- 
trinsic luminosity density ratio (Lisoo/i-LycOint = 7 (Siana 
et al. 2007; VlOb). In Sect. 3.2, the IGM transmission of 
Inoue et al. (2008,2011) is convolved with the adopted U- 
band filter. We perform calculations for three U-band fil- 
ters: F336W {F336W), u-band Stromgren (U-strom) and U- 
special available at the Large Binocular Telescope (LBC-U) 
(in general the results do not change if other filters with 
similar widths and wavelength coverages are considered). It 
is clear from Eq. (0 that high magnitude contrasts between 
the LyC and the 1500A rest-frame provide strong constraints 
on / e3C . We can rearrange Eq. [T] to give an estimate of the 
expected magnified flux at wavelengths smaller than the Ly- 
man limit as a function of / eS c and magnification factor /1 
(i.e., F^q 8 = /ii*Lyc): 
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As discussed in Inoue & Iwata (2008) and VlOb, the stochas- 
tic nature of the IGM absorption introduces large uncertain- 
ties in the estimate of f esc for a single line of sight. Stacking 
many galaxies among different lines of sight provides strong 
limits and reduces the variance due to IGM attenuation. 
Indeed, a deep U-band survey can provide firm constraints 
for L L* galaxies down to few percent of f esc by stacking 
tens of objects. This has been performed by VlOb with deep 
VLT/VIMOS U-band imaging (« 30,45 at la) of z > 3.4 
LBGs (see also Boutsia et al. (2011)). 

For L « L* galaxies the situation becomes challeng- 
ing. One reason is that the knowledge of the redshift is nec- 
essary to fix the observed wavelength position of the Lyman 
limit. Even though a deep U-band survey (e.g, mag-U~ 30 
at 1-sigma) can still probe f esc down to 20% by stacking 
about thirty 0.1L* galaxies, the redshift confirmation is not 
practicable with current facilities, or it is feasible only for 
a sub-class of sources like the Lya - emitters. The lensing 
magnification greatly facilitates the redshift measurement 
also for faint galaxies. Additionally it allows the investiga- 
tion of /esc down to stronger limits, by the U-band stacking 
and/or individual LyC detections. In the following we per- 
form Monte-Carlo (MC) simulations to derive predictions for 
/esc as a function of the U-band depth, lensing magnification 
and lensing cross-section. In what follows we adopt r and R 
to indicate the r-band magnitude (~ 1500 A rest-frame) for 
non-magnified and magnified sources, respectively. 



MAGNIFIED GALAXIES: 
SIMULATIONS 



MC 



We want now to estimate the feasibility of probing / eS c down 
to < 20% for L < 0.1L* LBGs. To this aim, from Eq. [2] it 
is possible to calculate as a function of f esc and the U-band 
depth the number of sources detectable in the LyC above a 
certain threshold and the signal to noise ratio in their LyC 
stack. We therefore perform MC simulations similar to those 
described in VlOb by assuming reasonable distributions for 
the quantities involved in Eq. [2] 
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Figure 1. MC simulations for magnified LBGs at 3 < z < 3.5 
with intrinsic L < 0.14L* (r > 26.5) and adopting a lensing 
cross-section of 10 arcmin 2 (see text for details). Top, middle and 
bottom panels show the results of our MC simulations for U-band 
depths of 29, 30 and 31 (1-sigma), respectively. Left panels: the 
expected S /N ratio in the LyC of the stacked fluxes (median and 
central 68% interval, solid and dashed lines, respectively). Filled 
circles are the S/N ratio (median and central 68% interval) in 
the "On/Off case" with 20% of the galaxies with f eS c =1 and 
80% with f esc =0. Horizontal dotted lines mark S/N=3. Right 
panels: number of direct detections in the LyC at S/N > 1 
(black solid line) and S/N > 2 (red dashed line). The dotted 
black lines report the number of galaxies in the faint magnitude 
bins from which they migrate to brighter apparent magnitudes 
(R < 26). The black solid line is the sum of the dotted lines. 
From top to bottom the dotted lines refer to the magnitude bins 
26.5 < r < 27.0, 27.0 < r < 27.5, etc. The red stars mark the 
expected number of direct detections at S/N > 2 in the "On/Off 
case" . 

3.1 The expected number of magnified LBGs 

The probability that a source galaxy is magnified by more 
than (j, can be expressed in terms of the probability den- 
sity P(u): P(> fi) = f°° P(fi)dn where P(/i) = -dP(> 
[i)/dfi. An interesting property of the lensing probability 
in the source plane is that P(> fx) tx /i~ 2 and therefore 
P(fi) oc fjT 3 for ^ » 1 as can be shown in particular cases 
and argued to be true in general (see Fig. 9 of Lima et 
al. 2010; Schneider et al. 1992). A maximum magnification 
(fimax) is imposed by the size of the source galaxies that 
we assume to have an intrinsic half light radius not smaller 
than 0.2kpc (assuming circular shape), i.e., u max = 400. 

Each galaxy cluster produces a region within which the 
magnification in the source plane is larger than fi m in- This 
area is an effective cross-section for lensing statistics. For 
our purposes, it is enough to perform a simple calculation 
for the cross section aiens'- 

Nelust 

ai ens ~ — — \_. {arcmin 2 ) (3) 



where 8e,i is the Einstein radius of a given cluster. Lenses 
with the largest Be are the best "cosmic telescopes", in 
particular we consider in our simulations the case with 
ciens = 10 arcmin 2 and fMmin = 1, corresponding approxi- 
mately to 4-5 massive galaxy clusters (see Sect. 4), the re- 
sults can be linearly rescaled to other areas. 

Adopting the z ~ 3 luminosity function parameters 
of Reddy & Steidel (2009) (a = -1.73, M* AB = -20.97, 
0* = 1.71 x 10" 3 Mpc -3 ), we can calculate the expected 
number of star-forming galaxies within the chosen area and 
redshift interval dz. A further important constraint is that 
the magnification must produce an apparent magnitude R 
brighter than the spectroscopic limit for redshift measure- 
ments (R < m spec ), necessary to fix exactly the wavelength 
position of the Lyman limit. We adopt m spec = 26 that 
corresponds to a relatively deep spectroscopic survey (e.g., 
Vanzella et al. 2009). Together with the limit on/j(< 400) 
defined above, the chosen m spec implies a limit on the inte- 
gration of the luminosity function, that is r = 32.5 (i.e., 
sources with r > 32.5 would require [i > 400 to have 
R < 26). It turns out that within an area ai ens = 10 arcmin 2 
and redshift interval 3.0 — 3.5, ~ 2300 LBGs are expected 
in the magnitude range 26.5 < r < 32.5, 30 of which have 
R < m apec = 26 (we refer to this sample as N spec = 30). 
Clearly, the intrinsically faint galaxies are the rarer cases ac- 
cessible from ground-based spectroscopy, being sources that 
need large magnification jj,. We anticipate here that the con- 
straints we discuss below on f esc are mainly dictated by 
sources belonging to the magnitude bin 26.5 < r < 28 (i.e., 
L — 0.05 — 0.16L*, on average 24 of the 30 galaxies). 



3.2 MC simulations of L < 0.1L* LBGs 

We run MC simulations on the area ui ens — 10 arcmin 2 
by varying for each galaxy the parameters involved in Eq. [2] 
and derive the ionised flux Fl/yQ S as a function of f eS c ■ The 
flux is then compared with the depth of the U-band survey 
under study. The procedure is described in details in VIOb. 
Briefly, a redshift is randomly extracted from the interval dz 
(uniformly) and associated to a galaxy with a given r mag- 
nitude extracted randomly from the magnitude range 26.5- 
32.5 (accordingly with the magnitude distribution provided 

IGM _ r IGM 

by the LF). The intergalactic transmission T LyC = e * 
has been derived from thousands of random realizations at 
the extracted redshift value and convolved with the adopted 
U-band filter. The lower limit z m i n of the redshift range dz 
is dictated by the adopted filter shape that probes LyC at 
z > 3, 3.2 and 3.3 for the three filters F336W, U-strom 
and LBC-U, respectively. The ideal upper limit would be 
Zmin + 0.1 such that the closest region blueward the Lyman 
limit is probed. However we relax to z m i„ + 0.5 by allowing 
more sources to be included (in order to reach N 3pec ~ 30 
in the three cases). The IGM prescription adopted here in- 
cludes all the intervening absorption systems and modu- 
lates properly the Fj/yc signal for increasing redshift. This 
is also supported by the recent direct detection of LyC at 
\rest < 830A (VIOb). 

The magnification factor (pi) is associated to each 
galaxy by extracting it randomly from the distribution 
P{n) ~ M~ 3 , as explained earlier in the text (we assume 
that galaxies are distributed uniformly over the sky). The 
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z „i„= 3 2 (U-strom) z mln = 3.0 (F336W) z mln =3.3 (LBC-U) 
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Figure 2. Signal to noise in the LyC of the stacking (top panels) 
and the fraction (%) of LyC detections with S/N higher than 1 
and 2 (bottom panels) are shown as a function of the U-band 
depth. Three different U-band filters have been used. In the three 
cases cr; ens = 10 arcmin 2 and N spec ~ 30. In all panels the three 
curves have been calculated for f esc =0.2, 0.15 and 0.10 from top 
to bottom, respectively. Horizontal dotted lines mark the one and 
three sigma levels. In the bottom panels solid and dashed lines 
refer to S/N > 1 and S/N > 2, respectively. The inset in the 
top panels show the adopted filter, with the dotted vertical line 
marking the position of the Lyman limit at the z m i n . 



with a U-band depth £ 30 (lc), the S/N of the stacking is 
> 1(3) for f esc > 0.05(0.18). The average number of indi- 
vidual detections (> 2a) is > 3 for f eS c > 0.10(0.30) if the 
U-band depth is 31(30). These detections are originated by 
magnification of sources with r > 26.5 that migrate from 
their original magnitude bins to the brighter lensed regime 
(R < 26). The contribution of each bin to the N spec sam- 
ple is shown as dotted lines in Figure U As a consistency 
check, and as expected, their number per magnitude bin 
flattens if the faint end slope of the LF is a = —3 (i.e., 
the growth of the number counts and the dimming of the 
lensing cross-section compensate each other). The flatten- 
ing of the number of direct detections as f esc increases (in 
the fixed f eac case, Figure [1} is due to IGM attenuation, i.e., 
the interception of Lyman limit systems (LLSs) or damped 
Lyman-alpha systems (DLAs) along the line of sight sup- 
presses the ionising flux even for high f esc values. 

Figure [5] shows the same quantities as a function of the 
U-band depth by fixing f esc to 10, 15 and 20% for all sources. 
It has been performed for three filters: F336W, U-strom and 
LBC-U. Clearly, the deeper U-band surveys give the best 
constraints to the faint galaxies, in particular those with 
mag-U> 30 are the more effective and increase rapidly as 
the mag-U= 31 limit is approached. The intermediate band 
filter like the u-stromgren is an optimal solution. It is worth 
noting that observing a massive cluster with 6e — 1 arcmin 
(such us MACS 0717.5+3745, see below), we expect ~ 10 
galaxies with R < 26 and r > 26.5. If all of them have 
fesc = 15%, then with a U-band (u-stromgren-like) depth 
of 30.5 we expect to detect 3 of them at S/N > 1 and 
a S/Nl v c — 3.0 for the stacking. Therefore, we expect to 
observe signal already with a single deep pointing (of ~60hr) 
with a 8-10m class telescope. 



fesc is explored in two different regimes of variability: (1) 
it is fixed to a constant value for all N sp ec sources span- 
ning the range 0.01-1.00, with steps of 0.01 (i.e., the average 

< fesc >= fesc of each galaxy) and (2) it is fixed to values 
or 1, assigned randomly such that the 80% of the N spe c 
sample have fesc =0 and 1 for the remaining sources, i.e., 

< fesc >= 0.2 with large variance (the "On/Off case"). 

In summary, the sample of galaxies derived from the LF 
in the considered volume (area of 10 arcmin 2 and dz= 0.5) 
and magnitude interval 26.5 < r < 32.5, has been extracted 
10000 times by varying / esc , e~ T * , (L 1500 / L LyC )int, 
(^i50o)obs and a (no dust has been considered). Figure [T] 
shows the results for the F336W case. The resulting median 
signal to noise ratio of the stacked LyC fluxes of the galax- 
ies with R < 26 is shown as a function of / esc and U-band 
depth (left panels). The signal to noise S/N L t y a g k is calcu- 
lated as E^T" KyrT 3 / Nspec] /[lo-u / V ' N sp ec] , where \o v 
is the depth of the U-band survey and N ap ec — 30, i.e., the 
number of surviving LBGs with R < 26. The number of di- 
rect LyC detections above the 2-sigma limit {Fl;g a > 2ou) 
is also reported in the same Figure Q] (right panels) . While 
the median S/NfjyQ k of the two / esc distributions, constant 
and "On/Off", are compatible if < fesc >= 0.2 (as expected 
from the S/Nx%§ k calculation), the variance is larger in the 
"On/Off" case. Also the number of direct detections are 
different between the two, because by definition we fix a 
maximum fraction of detectable sources. It turns out that 



4 THE CURRENT OBSERVATIONS 

As mentioned above, the interception of LLSs or DLAs 
makes IGM attenuation quite stochastic and severe for ion- 
ising radiation at z ^ 3, therefore the observations of mag- 
nified galaxies along different lines of sight are necessary to 
attenuate its effect. Assuming these absorbers are associated 
with individual structures with physical size of 50 kpc 
(physical) (e.g., Yajima et al. 2012), different lines of sight 
separated by more than ~ 10 arcsec will have a relatively 
small transverse correlation among high column density ab- 
sorbers, even along a single cluster. 

Currently, there are in the literature several massive 
clusters of galaxies characterised by relatively large lens- 
ing cross-sections, e.g., Be ~ 0.7 — 0.9 arcmin, like MACS 
0717.5+3745 (Zitrin et al. 2009), 1E0657-56 ("Bullet Clus- 
ter", Bradac et al. 2009) and A1689 (Broadhurst et al. 2005). 
Zitrin et al. (2012) recently measured the Einstein radius 
of 10000 galaxy clusters from the SDSS. In particular they 
identified ~40 candidates with 9e > 0.7 arcmin. The largest 
Einstein radius was 8e — 1-1 arcmin for the most massive 
cluster (for a source at z s = 2). The MACS survey (MAs- 
sive Cluster Survey) also provides a statistically complete 
sample of very X-ray luminous distant clusters of galaxies 
(0.3 < z < 0.7, Ebeling et al. 2010). Therefore, given the 
existing observations, the area of er; ens = 10 arcmin 2 con- 
sidered in this letter can be easily achieved with 5 clusters. 
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It is worth mentioning the ongoing Cluster Lensing And 
Supernova survey with Hubble (CLASH), a multi-treasury 
program that, in addition to other available ground-based 
observations is observing (524 HST orbits) 25 X-ray selected 
massive galaxy clusters with new HST panchromatic imag- 
ing capabilities (16 filters): the ACS and both the UVIS 
and IR channels of the WFC3. While 20 clusters were cho- 
sen to be X-ray selected, relatively relaxed clusters with no 
lensing-selection bias, the 5 additional were chosen to be 
high-magnification clusters with 6e > 0.6 arcmin (Postman 
et al. 2012). Even though the survey was not designed to pro- 
vide deep limits on LyC for faint galaxies, the homogeneous 
panchromatic photometry on a statistically significant sam- 
ple of massive clusters, together with magnification maps 
and high S/N spectroscopy on a number of highly magni- 
fied galaxies at z Si 3, will yield a very valuable testbed for 
probing LyC emission at the faint end of the galaxy LF. 



5 CONCLUSIONS 

We have explored the possibilities offered by gravitational 
lensing in constraining the ionizing emission from faint 
galaxies at z > 3, which has a great impact on studies of the 
reionization even at significantly higher redshifts (z > 6). 
We have shown that with current ground-based facilities and 
the magnification provided by known clusters of galaxies it 
is possible to constrain the / esc quantity down to 10-15% for 
faint luminosities (L £ 0.1L*). In particular, strong lensing 
provides the opportunity to make a significant step forward 
in this field in the pre-ELT era: (1) The magnification al- 
lows us to measure redshifts for very faint galaxies, i.e., to fix 
the wavelength position of the Lyman limit and put strong 
constraints on the LyC with U-band imaging for galaxies 
down to L ~ 0.05 — 0.1L*. (2) The magnification allows us 
to perform accurate spatial analysis of the LyC emission (if 
detected), down to hundreds of parsecs. If, on the one hand, 
the magnified area increases the probability of intercepting 
a foreground lower-z source that might mimic LyC emission 
(Vanzella et al. 2010a), on the other hand, the possible pres- 
ence of multiple images of the same background source will 
help to solve the problem, since a real LyC detection will be 
present in all the counterparts (if p are similar). The trans- 
verse separation of the light paths decreases rapidly as the 
source redshift is approached. Therefore the IGM attenua- 
tion is practically the same for all the multiple images. (3) 
It has been suggested that the f esc increases with decreasing 
luminosity (e.g., Yajima et al. 2011) and the ionizing back- 
ground is mainly produced by a large number of sub— L* 
galaxies. A nice feature of the present approach is to put a 
limit to the process of attributing the ionizing contribution 
to lower and lower luminosities with the increasing obser- 
vational depth of null results. This in fact would imply a 
steepening of the faint end of the LF that in principle can 
be excluded by the increasing number of highly magnified 
sources expected to be detected through lensing. 

We have shown that deep U-band observations (as in 
VlOb) of five massive galaxy clusters provide a sample of 
intrinsically faint (L £ 0.1L*) and magnified z ~ 3 galaxies 
(about 30) useful to investigate f eac down to 10-20%. An in- 
creasing number of well studied clusters for this study, with 



multi-band HST photometry, spectroscopy and lens model- 
ing, are becoming available as part of the CLASH project. 
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